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1. Introduction 



Relativistic mean-field theory seems to be a powerful tool in describing various aspects 
of nuclear structure [0, ^ 0. This theory provides an elegant and economical framework, 
in which properties of nuclear matter [||, ||], finite nuclei and neutron star, as well as 
the dynamics of heavy-ion collisions, can be calculated. Compared to conventional 
nonrelativistic approaches, relativistic models explicitly include mesonic degrees of 
freedom and describe nucleons as Dirac quasiparticles. Moreover, the spin-orbit 
interaction arises naturally from the Dirac-Lorenz structure of the effective Lagrangian. 
Nucleons interact in a relativistic covariant manner through the exchange of virtual 
mesons: the isoscalar scalar cr-meson, the isoscalar vector cu-meson, and the isovector 
vector p-meson. The model is based on the one-boson exchange description of the 
nucleon-nucleon interaction. One of the most interesting problems in intermediate 
energy physics which has been investigated with many approaches is an effective meson 
mass in nuclear medium. 



2. EfFective vector meson masses in nuclear medium. 

In this paper the properties of nuclear matter in high-density and neutron-rich regime 
are considered with the use of the relativistic mean-field approximation. This approach 
implies the interaction of nucleons through the exchange of meson fields, so the model 
considered here comprises: nucleons and scalar, vector-isoscalar and vector-isovector 
mesons. Consequently contributions coming from these components influence the 
properties of nucleon. In this paper the nonlinear vector-isoscalar self-interaction is 
dealt. Modification of this type was proposed by Bodmer in order to achieve good 
agreement with the Dirac-Briickner @ calculations at high densities. The first version of 
the p meson field introduction is of a minimal type without any nonlinearity and consists 
only of the coupling of this field to nucleons. This case is enlarged by the nonlinear 
vector-isoscalar and vector-isovector interaction and the vector-scalar interactions which 
modify the density dependence of the p mean field and the energy symmetry. Such an 
extension was inspired by the paper in which the authors indicate the existence of a 
relation between the neutron-rich skin of a heavy nucleus and the properties of a neutron 
star crust. In the quark meson coupling model [P] nucleon properties are modified by the 
meson coupling directly to the quarks and not to the nucleons ( paper by Li et.al. §). 
The RMF theory is an effective one, as coupling constants are determined by the bulk 
properties of nuclear matter such as saturation density, binding energy, compressibility 
and the symmetry energy. 

The starting point in this paper is the construction of the effective Lagrangian [l^, |ll], |l^ 



intended for application to the system described above. The scalar a, isoscalar-vector 
uj and isovector- vector p mesons are denoted by ip, cj^, 6^, respectively. The Lagrange 
density function for this model has the following form 
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- ^F^uF'"' + f^i^D^^ - ^{M - gN.v)i^. (1) 
The field tensors i?^^, VL^^^ F^^, and the covariant derivative are given by 

R1, = dX-d^bl + 9pea,chlhl, (2) 

Vt^u = d^Uu - d^Uf,, (3) 
= d^A, - d,A^, (4) 

D^ = d^ + ]^igphlT'' + ig^u^ + iQ^A^. (5) 

The potential function U{ip) possesses a polynomial form introduced by Boguta and 
Bodmer ||13| in order to get a correct value of the compressibility K of nuclear matter 
at saturation density 

u{v) = + \g2^'' + \g-iv^- (6) 

The nucleon field if) has a form of a column matrix composed of proton and neutron 
fields, respectively 

The nucleon mass is denoted by M whereas m^, M^, Mp are masses assigned to the 
meson fields. The parameters entering the Lagrangian function (|l|) are the coupling 
constants C3, g^, gp, g^a^ (Table ||) for meson fields, e^/in = 1/137 for photon field and 
the self-interacting coupling constants g2 and g^. The coupling constants g^^a and gpa in 
the QMC model [1^ are related to the constant g^a by the relations 
2 2 

guja — -gNa, gpa — -gNa- (8) 

The Lagrangian function (P includes also nonlinear terms 

\c,{uj,u;n' + A.{gpg^)\uj,uj^){b;b'^n (9) 
- g^M^^u.oo'' + ^gWuj.uj'^ - gpM.^b^b'^^ + ^g^^Xb^^ 

which affect remarkably properties of nucleons. The first term was proposed by Bodmer 
in order to achieve good agreement with the Dirac-Briickner calculations at high 
densities. All terms which involve the field change the density dependence of the 
symmetry energy and cause the change in the density dependence of the p field. The 
modification of the density dependence of the p mean field (|5|) can lead to the change of 
the neutron-skin thickness |^. This effect is observed in the Parity Radius Experiment 
(PREX) at the Jefferson Laboratory which aims to measure the neutron radius in ^°^Pb 
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via parity violating electron scattering [pTSl - 

The field equations for meson fields $^ = u;^, 6^} are 

+ {Mp-g^,^)gp,hlh^r 
- d^n^" = (M^ - g^^ipyoo" + csiuj^LO^P'' (11) 



+ 2A,{g^g^y{blb'^^)uj^-g^J: 



- D.R^-'^ = {Mp - g.^^fb'"' + 2A,ig^gpY{co^uj>^)b'''' - (?,(5,,3 ^3 • (12) 
Sources that appear in these equations are the baryon current 

Jb = ^7"^ (13) 
and the isospin current which exists only in the asymmetric matter 



J3 = 7;^rr'^ (14) 



1 
2 

The Maxwell equations are of the form 

-d^F"" = -e r 
with the electric current given by 

r = ^i^YQei^. (15) 

The timelike components of the baryon and isospin currants are the baryon density 
Pb = Jb the isospin density pa = The conserved baryon and isospin charges 
are given by the following relations: 

Qb= d^xip^ip = / d^x{i!:^ilJn + i'p'ipp), 



The electric charge is defined by 

Q = j d'x^+Q.ij = j d^xip;% , = i(/ - r^) = 1^ ° ° 

The Dirac equations for nucleons that are obtained from the Lagrangian function have 
the following form 

i^Dpij - (M - g^Mi^ = 0. (16) 
The physical system is defined by the thermodynamic potential ||T6| 

Q = -A;T/nTr(e-^(^-^^^-^^^^)) (17) 
where H stands for the Hamiltonian and is defined as 

H = J2 [ d^x{do<^Arf^ - C}. (18) 

A J 
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71"^ = are momenta connected with fields $a = {v^, which denote 

all fields considered in this model. The averaged charges can be obtained from the 
thermodynamic potential with the use of the relations 

In this paper the variational method based on the Feynman-Bogolubov inequality [|r 



is incorporated (see more details in p8| ) 

n<ni = ^]o("^e//)+ < h- Ho>o . (20) 

Qq is the thermodynamic potential of the trial system which as effectively free 
quasiparticle system is described by the Lagrange function 

--^:.G'''' + \Y.<f..W (21) 

a 

+ ip{i-f^D^ - meff^N)^^ 
with the field tensor now defined by 

= dX - ^X 

and 

The fields $a can be written as a sum of two components, the effectively free 
quasiparticle fields = Wi^iiX ^'^^ classical boson condensates = cr,w,b 

^A=^A + U, (22) 
which in the case of particular meson fields results in 

ip = Tp + a (23) 

ujf, = uJ^, + Wf,, = 6^flW (24) 

P; = ? + /5^, /?^ = 5»'35^,o& (25) 

Both the boson and fermion masses {mB,mF) as well as the C,a = {c, w, b} fields are 
treated as the variational parameters in the effective potential. The covariant derivative 
for the trial system is 

D^ = d^ + ^igpPy + ig^wp (26) 

This introduces the fermion interaction with homogenous boson condensate w^, The 
fermion quasiparticle obeys the Dirac equation 

(27'^D^-me//,;v)^ = (27) 
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Figure 1. 




The masses entering the Lagrangian function LQ{meff) of the trial system calculated 
from the extremum conditions 

dQi dQi dQi dQi 
— = — = — = — = 



^^efS^a ^^eSUu ^^eff,p dm^ff^N 



now are the effective once. They are given by the following equations 

"^e//,<x = ml + 3^3 < >0 +262O- + 363^2 (28) 

- < ^t^^^ >o) - g%ip^+ < ^5"" >o), 

Kff,u^ = - Qau^fyf + gL < >0 +C3(W^+ < UJ^U)^ >o) (29) 

+ 2K{g,g^)\\?^ < >o) 
mlfj^p^, = {Mp - gp^af + gl„ < if >o +2{gpg^fK{w^+ < IZpUj^ >o) 
-^>^(l-M-^^E<5^">o (30) 

meff,N = M - gN„a. (31) 

In the presence of the p field condensation h the p mass nieff^p^a is splitting. This is a 
result of the 'nonabelian' character of the SU{2) symmetry which is explicitely broken 
by the Mp term. In the homogeneous case the free energy reaches the minimum at a. 
The nucleon effective mass meff,N dependence on the baryon density is presented 
in Fig. |1|. Results are shown for the TMl and NL3 parameter sets without additional 
coupling and for NL3 supplemented by the vector-scalar coupling. 

The meson effective masses as functions of the baryon density hb are presented in 
Figs. H and This figure provides a comparison of the meson effective masses obtained 
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Figure 2. The mesons uj and p effective mass as the function of the baryon density 
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Figure 3. The mesons lu and a effective mass as the function of the baryon density 
nB {fm~^). 

for different values of the coupling constants. Straight dotted lines show the u and p 
masses in the simplest model with no additional nonlinear terms. When this model 
is enlarged by the nonlinear vector-isoscalar and vector-isovector and the vector-scalar 
interactions the effective masses are obtained. As baryon number increases the u meson 
effective mass increases and then become smaller. The p meson effective mass increases 
with the increasing baryon number. 

The symmetry energy density (for T = 0) Eg = Es/ub can be evaluated taking into 
account the contributions of the p meson field to the energy of the system 

= l^lff^ = liK- 9pa(yf + A,{gpg^fw%^ (32) 
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Using the equation of motion the static homogeneous solution for the p field takes the 
following form 

b = -4^Q. = (33) 

S = {rin- np)/nB (34) 

and is presented in Fig. H. Taking into account the form of the equations (|3^) and (^) 
the symmetry energy now is given by 
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Parameter 


NL3 


TMl 


QMC 


Eo (MeV) 


-15.92 


-16.26 


-16.65 


kF,o ifm-^) 


1.24 


1.29 


1.297 


Po ifm-^) 


0.14 


0.145 


0.147 


K (MeV) 


170.24 


281.53 


284.71 



Table 1 . The nucleon symmetric phase properties (the binding energy Eq , the Fermi 
momentum, the density and the compresibiHty K) at the saturation point. 



In the case of no additional coupling constants the result reproduces that of the simplest 
model without any nonlinearities. Including the u — p and the vector-scalar coupling the 
Es is no longer linear with baryon density n^. Including the u — p and the vector-scalar 
coupling the Es is no longer linear with baryon number density Ub- For the purpose 



of the present study, we choose two basic RMF parameterizations (Table y) : NL3 ||22 
TMl and one based on QMC approach ||l4|. All parameters are chosen in such a 
form to reproduce the nucleon symmetric phase properties (the binding energy Eq, the 
Fermi momentum, the density and the compresibility K) at the satutation point (Table 

i- 

The case of static nonhomogeneous nuclear matter configuration can be described 
by the Euler-Lagrange equations (|T0| , 0, |ll], |l^) taking the quantum average of these 



equations with respect to the effective quasiparticle system Hq which in turn is obtained 

with the use of the effective quasiparticle Lagrangian £o("^e//)- The results for a 
spherically symmetric system become 

Aa(r) = (35) 

Aw{r) =mlff^^w{r) - g^Qe, (36) 

Abir)=mlff^pbir)-gpQ3 (37) 

where the effective potential for the scalar meson is given by 

Ueff = {g2 < Tp^ >o +^^aM^(w^+ < uj^uj^ >o) + gp^Mp{b^+ < >o))a 



+ ^{ml + 3(73 < >o -9L{w^+ < ^M^^ >o) - glAb^+ < >o)) c^' 



+ l92(r^ + ^gsa\ (38) 
The construction of Ueff is such that 

2 _ d^Ueff 
^eff,. - -g^- 

The p condensate has a mass (rrieff^p = rrieff^p^s at T = 0) 

mlfjp = {Mp - gp^a f + gl^ < >q +2{gpg^YA^{w^+ < UpU" >o) 
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Figure 6. Profile of the ui meson field inside the ^osp^ nucleus {g^a = 9 pa = case). 



In the obtained results the thermal contributions from meson fields are present, however, 
they are not involved in calculations since meson masses are too large to be relevant. 
The Dirac equation for nucleons (i=p,n) has the following form 

{-zaV + V{r) + (3[M - 5(r)]}^, = e.^j, (40) 

with the potentials: 

V{r) = g^w{r) + Qp a%{r) + eQeAQ{r), 

Sir) = gNaCrir). 

In the non-honomegenous nuclear matter distribution the effective masses are space 
dependent and act rather like external potentials. The self-consistent equations 
|36| , |37| , |40|) , namely the Dirac equation with potential terms for the nucleons and the 
Klein-Gordon type equations with sources for the mesons and the photon are the base 
for the numerical description of nuclei (l^. The simplified solutions (cf. Fig. |^ and 
Fig.U) can be obtained assuming the Fermi shape for the neutron and proton density 
distribution 0. Spherically symmetric solutions for the u and p meson fields for the 
^°^Pb nucleus {g^^(J = gp„ = case) are presented in Fig. ^ and Fig. ^ More realistic 
calculations, which will be the subject of future investigations, can be obtained with the 



use of the self-consistent Thomas-Fermi approximation |^, |^. In the case of the uj 
meson the influence of the u — p coupling is rather small whereas there exists a marked 
difference in the p meson field profile when the u — p coupling is present (cf. (Fig. |^)). 



3. Conclusions 



The model presented here includes the original Walecka model with the nucleon fields 
and the kinetic and mass terms for scalar and vector mesons. Scalar nonlinear terms 
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r( rm ) 



Figure 7. Profile of the p meson field inside the ^"^Pb nucleus {guia — 9 pa — case). 



9pM 


A^[0] 




{Rn-Rp) (/m)0 


me//,^ {MeV) 


^eff,p {MeV) 


8.92 





5.460 


0.280 


782.5 


763.0 


10.75 


0.025 


5.469 


0.209 


785.4 


1141.6 



Table 2. The meson effective masses inside the ^"^Pb nucleus {guia — 9 pa — case). 



were proposed by Boguta and Bodmer. Boguta and Price introduced quartic vector self- 
interaction. The chosen form of the Lagrange function which was supplemented by the 
additional omega-rho coupling was inspired by the work by Horowitz and Piekarewicz 
and by the QMC model [l^. This additional nonlinear meson interaction can be 
explained in terms of effective meson masses that are modified in the nuclear medium. 
The additional nonlinear couplings Ay and Quia and Qp^ change the density dependence 
of the omega and rho fields, whereas Ay itself influences the density dependence of 
the symmetry energy. Taking into account these additional interaction terms the 
selfconsistent equations for vector mesons were obtained. The nucleon and scalar self- 
interaction generates different from zero expectation value of the scalar fleld and the 
effective nucleon mass rrieff^N- The mutual interactions of vector meson flelds give as 
a result the effective meson masses meff^uj and rrieff^p. The appearance of the effective 
meson masses changes the energy symmetry of the system and has an influence on the 
behaviour of the p fleld. All calculations were performed with the use of NL3 |^ and 
TMl (2^ parameter sets. They not only correctly reproduce nuclear properties but 
also are consistent with the idea of naturalness as well. The force of NL3 stems from 
the flt including exotic nuclei, neutron radii, and information on giant resonances. The 
parameter set named TMl has been applied to deformed nuclei, triaxal deformations, 
nuclear excitations fl^, |2j, |2^. The last column mark by QMC represents parameters 



produced by the Quark Meson Coupling model (QMC) |T4| . 
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Parameter 


NL3 22 


NL3 


NL3 


TMl 2^ 


OMC [Il3l 


M 


938 MeV 


938 MeV 


938 MeV 


938 MeV" 


939 MeV 




795.4 Me\^ 


795.4 MeV 


795.4 MeV 


783 MeV 


783 MeV 


±y± p 


763 MeV 


763 MeV 


763 MeV 


770 MeV 


770 MeV" 


1 1 ug 


492 Mel^ 


492 MeV 


492 Me\/ 


511.2 Mel^ 


550 MeV 


Qo = /t:/2 


12.17 /m"^ 


12.17 fm^'^ 


12.17/m"^ 


7.23 /m"^ 


9.62 /m"^ 


Gfo = A/6 


-36.259 


-36.259 


-36.259 


0.6183 


61.24 


QAJ^ — Qa 


10.138 


10.138 


10.138 


10.0289 


6.9667 


Q< ! 


13.285 


13.285 


13.285 


12.6139 


5.8561 




8.941 


9.214 


10.751 


9.2644 


8.3803 













71.3075 










0.005 


0.025 















2/3 gNcT 





0.6617 


9pa 







2/3 gNcT 





0.6617 



Table 3. Parameter sets for the Lagrangian (^. 
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